HIHISIV: a database of gene expression in HIV and SIV host immune response.

In the battle of the host against lentiviral pathogenesis, the immune response is crucial. However, several questions remain unanswered about the interaction with different viruses and their influence on disease progression. The simian immunodeficiency virus (SIV) infecting nonhuman primates (NHP) is widely used as a model for the study of the human immunodeficiency virus (HIV) both because they are evolutionarily linked and because they share physiological and anatomical similarities that are largely explored to understand the disease progression. The HIHISIV database was developed to support researchers to integrate and evaluate the large number of transcriptional data associated with the presence/absence of the pathogen (SIV or HIV) and the host response (NHP and human). The datasets are composed of microarray and RNA-Seq gene expression data that were selected, curated, analyzed, enriched, and stored in a relational database. Six query templates comprise the main data analysis functions and the resulting information can be downloaded. The HIHISIV database, available at  https://hihisiv.github.io , provides accurate resources for browsing and visualizing results and for more robust analyses of pre-existing data in transcriptome repositories.

Measuring criticality in control of complex biological networks.

Recent controllability analyses have demonstrated that driver nodes tend to be associated to genes related to important biological functions as well as human diseases. While researchers have focused on identifying critical nodes, intermittent nodes have received much less attention. Here, we propose a new efficient algorithm based on the Hamming distance for computing the importance of intermittent nodes using a Minimum Dominating Set (MDS)-based control model. We refer to this metric as criticality. The application of the proposed algorithm to compute criticality under the MDS control framework allows us to unveil the biological importance and roles of the intermittent nodes in different network systems, from cellular level such as signaling pathways and cell-cell interactions such as cytokine networks, to the complete nervous system of the nematode worm C. elegans. Taken together, the developed computational tools may open new avenues for investigating the role of intermittent nodes in many biological systems of interest in the context of network control.

A genome-scale metabolic model of parasitic whipworm.

Genome-scale metabolic models are widely used to enhance our understanding of metabolic features of organisms, host-pathogen interactions and to identify therapeutics for diseases. Here we present iTMU798, the genome-scale metabolic model of the mouse whipworm Trichuris muris. The model demonstrates the metabolic features of T. muris and allows the prediction of metabolic steps essential for its survival. Specifically, that Thioredoxin Reductase (TrxR) enzyme is essential, a prediction we validate in vitro with the drug auranofin. Furthermore, our observation that the T. muris genome lacks gsr-1 encoding Glutathione Reductase (GR) but has GR activity that can be inhibited by auranofin indicates a mechanism for the reduction of glutathione by the TrxR enzyme in T. muris. In addition, iTMU798 predicts seven essential amino acids that cannot be synthesised by T. muris, a prediction we validate for the amino acid tryptophan. Overall, iTMU798 is as a powerful tool to study not only the T. muris metabolism but also other Trichuris spp. in understanding host parasite interactions and the rationale design of new intervention strategies.

Special Issue: New Advances in Bioinformatics and Biomedical Engineering Using Machine Learning Techniques, IWBBIO-2022.

Bioinformatics is revolutionizing Biomedicine in the way we treat and diagnose pathologies related to biological manifestations resulting from variations or mutations of our DNA [...].

Functional selectivity of Receptor Tyrosine Kinases regulates distinct cellular outputs.

Functional selectivity refers to the activation of differential signalling and cellular outputs downstream of the same membrane-bound receptor when activated by two or more different ligands. Functional selectivity has been described and extensively studied for G-protein Coupled Receptors (GPCRs), leading to specific therapeutic options for dysregulated GPCRs functions. However, studies regarding the functional selectivity of Receptor Tyrosine Kinases (RTKs) remain sparse. Here, we will summarize recent data about RTK functional selectivity focusing on how the nature and the amount of RTK ligands and the crosstalk of RTKs with other membrane proteins regulate the specificity of RTK signalling. In addition, we will discuss how structural changes in RTKs upon ligand binding affects selective signalling pathways. Much remains to be known about the integration of different signals affecting RTK signalling specificity to orchestrate long-term cellular outcomes. Recent advancements in omics, specifically quantitative phosphoproteomics, and in systems biology methods to study, model and integrate different types of large-scale omics data have increased our ability to compare several signals affecting RTK functional selectivity in a global, system-wide fashion. We will discuss how such methods facilitate the exploration of important signalling hubs and enable data-driven predictions aiming at improving the efficacy of therapeutics for diseases like cancer, where redundant RTK signalling pathways often compromise treatment efficacy.

Fluxomics reveals cellular and molecular basis of increased renal ammoniagenesis.

The kidney plays a critical role in excreting ammonia during metabolic acidosis and liver failure. The mechanisms behind this process have been poorly explored. The present study combines results of in vivo experiments of increased total ammoniagenesis with systems biology modeling, in which eight rats were fed an amino acid-rich diet (HD group) and eight a normal chow diet (AL group). We developed a method based on elementary mode analysis to study changes in amino acid flux occurring across the kidney in increased ammoniagenesis. Elementary modes represent minimal feasible metabolic paths in steady state. The model was used to predict amino acid fluxes in healthy and pre-hyperammonemic conditions, which were compared to experimental fluxes in rats. First, we found that total renal ammoniagenesis increased from 264 ± 68 to 612 ± 87 nmol (100 g body weight)-1 min-1 in the HD group (P = 0.021) and a concomitated upregulation of NKCC2 ammonia and other transporters in the kidney. In the kidney metabolic model, the best predictions were obtained with ammonia transport as an objective. Other objectives resulting in a fair correlation with the measured fluxes (correlation coefficient >0.5) were growth, protein uptake, urea excretion, and lysine and phenylalanine transport. These predictions were improved when specific gene expression data were considered in HD conditions, suggesting a role for the mitochondrial glycine pathway. Further studies are needed to determine if regulation through the mitochondrial glycine pathway and ammonia transporters can be modulated and how to use the kidney as a therapeutic target in hyperammonemia.

Spatially resolved phosphoproteomics reveals fibroblast growth factor receptor recycling-driven regulation of autophagy and survival.

Receptor Tyrosine Kinase (RTK) endocytosis-dependent signalling drives cell proliferation and motility during development and adult homeostasis, but is dysregulated in diseases, including cancer. The recruitment of RTK signalling partners during endocytosis, specifically during recycling to the plasma membrane, is still unknown. Focusing on Fibroblast Growth Factor Receptor 2b (FGFR2b) recycling, we reveal FGFR signalling partners proximal to recycling endosomes by developing a Spatially Resolved Phosphoproteomics (SRP) approach based on APEX2-driven biotinylation followed by phosphorylated peptides enrichment. Combining this with traditional phosphoproteomics, bioinformatics, and targeted assays, we uncover that FGFR2b stimulated by its recycling ligand FGF10 activates mTOR-dependent signalling and ULK1 at the recycling endosomes, leading to autophagy suppression and cell survival. This adds to the growing importance of RTK recycling in orchestrating cell fate and suggests a therapeutically targetable vulnerability in ligand-responsive cancer cells. Integrating SRP with other systems biology approaches provides a powerful tool to spatially resolve cellular signalling.

A calpain-6/YAP axis in sarcoma stem cells that drives the outgrowth of tumors and metastases.

Sarcomas include cancer stem cells, but how these cells contribute to local and metastatic relapse is largely unknown. We previously showed the pro-tumor functions of calpain-6 in sarcoma stem cells. Here, we use an osteosarcoma cell model, osteosarcoma tissues and transcriptomic data from human tumors to study gene patterns associated with calpain-6 expression or suppression. Calpain-6 modulates the expression of Hippo pathway genes and stabilizes the hippo effector YAP. It also modulates the vesicular trafficking of ß-catenin degradation complexes. Calpain-6 expression is associated with genes of the G2M phase of the cell cycle, supports G2M-related YAP activities and up-regulated genes controlling mitosis in sarcoma stem cells and tissues. In mouse models of bone sarcoma, most tumor cells expressed calpain-6 during the early steps of tumor out-growth. YAP inhibition prevented the neoformation of primary tumors and metastases but had no effect on already developed tumors. It could even accelerate lung metastasis associated with large bone tumors by affecting tumor-associated inflammation in the host tissues. Our results highlight a specific mechanism involving YAP transcriptional activity in cancer stem cells that is crucial during the early steps of tumor and metastasis outgrowth and that could be targeted to prevent sarcoma relapse.

Data integration and mechanistic modelling for breast cancer biology: Current state and future directions.

Breast cancer is one of the most common cancers threatening women worldwide. A limited number of available treatment options, frequent recurrence, and drug resistance exacerbate the prognosis of breast cancer patients. Thus, there is an urgent need for methods to investigate novel treatment options, while taking into account the vast molecular heterogeneity of breast cancer. Recent advances in molecular profiling technologies, including genomics, epigenomics, transcriptomics, proteomics and metabolomics data, enable approaching breast cancer biology at multiple levels of omics interaction networks. Systems biology approaches, including computational inference of 'big data' and mechanistic modelling of specific pathways, are emerging to identify potential novel combinations of breast cancer subtype signatures and more diverse targeted therapies.

SubcellulaRVis: a web-based tool to simplify and visualise subcellular compartment enrichment.

Cells contain intracellular compartments, including membrane-bound organelles and the nucleus, and are surrounded by a plasma membrane. Proteins are localised to one or more of these cellular compartments; the correct localisation of proteins is crucial for their correct processing and function. Moreover, proteins and the cellular processes they partake in are regulated by relocalisation in response to various cellular stimuli. High-throughput 'omics experiments result in a list of proteins or genes of interest; one way in which their functional role can be understood is through the knowledge of their subcellular localisation, as deduced through statistical enrichment for Gene Ontology Cellular Component (GOCC) annotations or similar. We have designed a bioinformatics tool, named SubcellulaRVis, that compellingly visualises the results of GOCC enrichment for quick interpretation of the localisation of a group of proteins (rather than single proteins). We demonstrate that SubcellulaRVis precisely describes the subcellular localisation of gene lists whose locations have been previously ascertained. SubcellulaRVis can be accessed via the web (http://phenome.manchester.ac.uk/subcellular/) or as a stand-alone app (https://github.com/JoWatson2011/subcellularvis). SubcellulaRVis will be useful for experimental biologists with limited bioinformatics expertise who want to analyse data related to protein (re)localisation and location-specific modules within the intracellular protein network.

Cytosolic fumarase acts as a metabolic fail-safe for both high and low temperature acclimation of Arabidopsis thaliana.

Plants acclimate their photosynthetic capacity (Pmax) in response to changing environmental conditions. In Arabidopsis thaliana, photosynthetic acclimation to cold requires the accumulation of the organic acid fumarate, catalysed by a cytosolically localized fumarase, FUM2. However, the role of this accumulation is currently unknown. Here, we use an integrated experimental and modelling approach to examine the role of FUM2 and fumarate across the physiological temperature range. We have studied three genotypes: Col-0; a fum2 mutant in a Col-0 background; and C24, an accession with reduced FUM2 expression. While low temperature causes an increase in Pmax in the Col-0 plants, this parameter decreases following exposure of plants to 30 °C for 7 d. Plants in which fumarate accumulation is partially (C24) or completely (fum2) abolished show a reduced acclimation of Pmax across the physiological temperature range (i.e. Pmax changes less in response to changing temperature). To understand the role of fumarate accumulation, we have adapted a reliability engineering technique, Failure Mode and Effect Analysis (FMEA), to formalize a rigorous approach for ranking metabolites according to the potential risk that they pose to the metabolic system. FMEA identifies fumarate as a low-risk metabolite, while its precursor, malate, is shown to be high risk and liable to cause system instability. We propose that the role of FUM2 is to provide a fail-safe in order to control malate concentration, maintaining system stability in a changing environment. We suggest that FMEA is a technique that is not only useful in understanding plant metabolism but can also be used to study reliability in other systems and synthetic pathways.

COVID-19 vaccination strategies depend on the underlying network of social interactions.

Since the onset of the coronavirus disease 2019 (COVID-19) pandemic, different mitigation and management strategies limiting economic and social activities have been implemented across many countries. Despite these strategies, the virus continues to spread and mutate. As a result, vaccinations are now administered to suppress the pandemic. Current COVID-19 epidemic models need to be expanded to account for the change in behaviour of new strains, such as an increased virulence and higher transmission rate. Furthermore, models need to account for an increasingly vaccinated population. We present a network model of COVID-19 transmission accounting for different immunity and vaccination scenarios. We conduct a parameter sensitivity analysis and find the average immunity length after an infection to be one of the most critical parameters that define the spread of the disease. Furthermore, we simulate different vaccination strategies and show that vaccinating highly connected individuals first is the quickest strategy for controlling the disease.

Mechanical loading activates the YAP/TAZ pathway and chemokine expression in the MLO-Y4 osteocyte-like cell line.

Osteocytes are mechanosensitive cells that control bone remodeling in response to mechanical loading. To date, specific signaling pathways modulated by mechanical loading in osteocytes are not well understood. Yes associated protein (YAP) and transcriptional coactivator with PDZ-binding motif (TAZ), the main effectors of the Hippo pathway, are reported to play a role in mechanotransduction and during osteoblastogenesis. Here, we hypothesized that YAP/TAZ signaling mediates osteocyte mechanosensing to target genes of the bone remodeling process. We aimed to investigate the contribution of YAP/TAZ in modulating the gene expression in an osteocyte-like cell line MLO-Y4. We developed a 3D osteocyte compression culture model from an MLO-Y4 osteocyte cell line embedded in concentrated collagen hydrogel. 3D-mechanical loading led to the increased expression of mechanosensitive genes and a subset of chemokines, including M-csf, Cxcl1, Cxcl2, Cxcl3, Cxcl9, and Cxcl10. The transcription regulators YAP and TAZ translocated to the nucleus and upregulated their target genes and proteins. RNAseq analysis revealed that YAP/TAZ knockdown mediated the regulation of several genes including osteocyte dendrite formation. Use of YAP/TAZ knockdown partially blunted the increase in M-csf and Cxcl3 levels in response to MLO-Y4 compression. These findings demonstrate that YAP/TAZ signaling is required for osteocyte-like cell mechano-transduction, regulates the gene expression profiles and controls chemokine expression.

Transcriptional Profile of the Industrial Hybrid Saccharomyces pastorianus Reveals Temperature-Dependent Allele Expression Bias and Preferential Orthologous Protein Assemblies.

Saccharomyces pastorianus is a natural yeast evolved from different hybridization events between the mesophilic S. cerevisiae and the cold-tolerant S. eubayanus. This complex aneuploid hybrid carries multiple copies of the parental alleles alongside specific hybrid genes and encodes for multiple protein isoforms which impart novel phenotypes, such as the strong ability to ferment at low temperature. These characteristics lead to agonistic competition for substrates and a plethora of biochemical activities, resulting in a unique cellular metabolism. Here, we investigated the transcriptional signature of the different orthologous alleles in S. pastorianus during temperature shifts. We identified temperature-dependent media-independent genes and showed that 35% has their regulation dependent on extracellular leucine uptake, suggesting an interplay between leucine metabolism and temperature response. The analysis of the expression of ortholog parental alleles unveiled that the majority of the genes expresses preferentially one parental allele over the other and that S. eubayanus-like alleles are significantly over-represented among the genes involved in the cold acclimatization. The presence of functionally redundant parental alleles may impact on the nature of protein complexes established in the hybrid, where both parental alleles are competing. Our expression data indicate that the majority of the protein complexes investigated in the hybrid are likely to be either exclusively chimeric or unispecific and that the redundancy is discouraged, a scenario that fits well with the gene balance hypothesis. This study offers the first overview of the transcriptional pattern of S. pastorianus and provides a rationalization for its unique industrial traits at the expression level.

Using Multilayer Heterogeneous Networks to Infer Functions of Phosphorylated Sites.

Mass spectrometry-based quantitative phosphoproteomics has become an essential approach in the study of cellular processes such as signaling. Commonly used methods to analyze phosphoproteomics datasets depend on generic, gene-centric annotations such as Gene Ontology terms, which do not account for the function of a protein in a particular phosphorylation state. Analysis of phosphoproteomics data is hampered by a lack of phosphorylated site-specific annotations. We propose a method that combines shotgun phosphoproteomics data, protein-protein interactions, and functional annotations into a heterogeneous multilayer network. Phosphorylation sites are associated to potential functions using a random walk on the heterogeneous network (RWHN) algorithm. We validated our approach against a model of the MAPK/ERK pathway and functional annotations from PhosphoSitePlus and were able to associate differentially regulated sites on the same proteins to their previously described specific functions. We further tested the algorithm on three previously published datasets and were able to reproduce their experimentally validated conclusions and to associate phosphorylation sites with known functions based on their regulatory patterns. Our approach provides a refinement of commonly used analysis methods and accurately predicts context-specific functions for sites with similar phosphorylation profiles.

A Holistic Approach to Study Photosynthetic Acclimation Responses of Plants to Fluctuating Light.

Plants in natural environments receive light through sunflecks, the duration and distribution of these being highly variable across the day. Consequently, plants need to adjust their photosynthetic processes to avoid photoinhibition and maximize yield. Changes in the composition of the photosynthetic apparatus in response to sustained changes in the environment are referred to as photosynthetic acclimation, a process that involves changes in protein content and composition. Considering this definition, acclimation differs from regulation, which involves processes that alter the activity of individual proteins over short-time periods, without changing the abundance of those proteins. The interconnection and overlapping of the short- and long-term photosynthetic responses, which can occur simultaneously or/and sequentially over time, make the study of long-term acclimation to fluctuating light in plants challenging. In this review we identify short-term responses of plants to fluctuating light that could act as sensors and signals for acclimation responses, with the aim of understanding how plants integrate environmental fluctuations over time and tailor their responses accordingly. Mathematical modeling has the potential to integrate physiological processes over different timescales and to help disentangle short-term regulatory responses from long-term acclimation responses. We review existing mathematical modeling techniques for studying photosynthetic responses to fluctuating light and propose new methods for addressing the topic from a holistic point of view.

Combinatorial mathematical modelling approaches to interrogate rear retraction dynamics in 3D cell migration.

Cell migration in 3D microenvironments is a complex process which depends on the coordinated activity of leading edge protrusive force and rear retraction in a push-pull mechanism. While the potentiation of protrusions has been widely studied, the precise signalling and mechanical events that lead to retraction of the cell rear are much less well understood, particularly in physiological 3D extra-cellular matrix (ECM). We previously discovered that rear retraction in fast moving cells is a highly dynamic process involving the precise spatiotemporal interplay of mechanosensing by caveolae and signalling through RhoA. To further interrogate the dynamics of rear retraction, we have adopted three distinct mathematical modelling approaches here based on (i) Boolean logic, (ii) deterministic kinetic ordinary differential equations (ODEs) and (iii) stochastic simulations. The aims of this multi-faceted approach are twofold: firstly to derive new biological insight into cell rear dynamics via generation of testable hypotheses and predictions; and secondly to compare and contrast the distinct modelling approaches when used to describe the same, relatively under-studied system. Overall, our modelling approaches complement each other, suggesting that such a multi-faceted approach is more informative than methods based on a single modelling technique to interrogate biological systems. Whilst Boolean logic was not able to fully recapitulate the complexity of rear retraction signalling, an ODE model could make plausible population level predictions. Stochastic simulations added a further level of complexity by accurately mimicking previous experimental findings and acting as a single cell simulator. Our approach highlighted the unanticipated role for CDK1 in rear retraction, a prediction we confirmed experimentally. Moreover, our models led to a novel prediction regarding the potential existence of a 'set point' in local stiffness gradients that promotes polarisation and rapid rear retraction.

Transgenic inhibition of interleukin-6 trans-signaling does not prevent skeletal pathologies in mucolipidosis type II mice.

Severe skeletal alterations are common symptoms in patients with mucolipidosis type II (MLII), a rare lysosomal storage disorder of childhood. We have previously reported that progressive bone loss in a mouse model for MLII is caused by an increased number of bone-resorbing osteoclasts, which is accompanied by elevated expression of the cytokine interleukin-6 (IL-6) in the bone microenvironment. In the present study we addressed the question, if pharmacological blockade of IL-6 can prevent the low bone mass phenotype of MLII mice. Since the cellular IL-6 response can be mediated by either the membrane-bound (classic signaling) or the soluble IL-6 receptor (trans-signaling), we first performed cell culture assays and found that both pathways can increase osteoclastogenesis. We then crossed MLII mice with transgenic mice expressing the recombinant soluble fusion protein sgp130Fc, which represents a natural inhibitor of IL-6 trans-signaling. By undecalcified histology and bone-specific histomorphometry we found that high circulating sgp130Fc levels do not affect skeletal growth or remodeling in wild-type mice. Most importantly, blockade of IL-6 trans-signaling did neither reduce osteoclastogenesis, nor increase bone mass in MLII mice. Therefore, our data clearly demonstrate that the bone phenotype of MLII mice cannot be corrected by blocking the IL-6 trans-signaling.

Metabolic flux from the chloroplast provides signals controlling photosynthetic acclimation to cold in Arabidopsis thaliana.

Photosynthesis is especially sensitive to environmental conditions, and the composition of the photosynthetic apparatus can be modulated in response to environmental change, a process termed photosynthetic acclimation. Previously, we identified a role for a cytosolic fumarase, FUM2 in acclimation to low temperature in Arabidopsis thaliana. Mutant lines lacking FUM2 were unable to acclimate their photosynthetic apparatus to cold. Here, using gas exchange measurements and metabolite assays of acclimating and non-acclimating plants, we show that acclimation to low temperature results in a change in the distribution of photosynthetically fixed carbon to different storage pools during the day. Proteomic analysis of wild-type Col-0 Arabidopsis and of a fum2 mutant, which was unable to acclimate to cold, indicates that extensive changes occurring in response to cold are affected in the mutant. Metabolic and proteomic data were used to parameterize metabolic models. Using an approach called flux sampling, we show how the relative export of triose phosphate and 3-phosphoglycerate provides a signal of the chloroplast redox state that could underlie photosynthetic acclimation to cold.

HybridMine: A Pipeline for Allele Inheritance and Gene Copy Number Prediction in Hybrid Genomes and Its Application to Industrial Yeasts.

Genome-scale computational approaches are opening opportunities to model and predict favorable combination of traits for strain development. However, mining the genome of complex hybrids is not currently an easy task, due to the high level of redundancy and presence of homologous. For example, Saccharomyces pastorianus is an allopolyploid sterile yeast hybrid used in brewing to produce lager-style beers. The development of new yeast strains with valuable industrial traits such as improved maltose utilization or balanced flavor profiles are now a major ambition and challenge in craft brewing and distilling industries. Moreover, no genome annotation for most of these industrial strains have been published. Here, we developed HybridMine, a new user-friendly, open-source tool for functional annotation of hybrid aneuploid genomes of any species by predicting parental alleles including paralogs. Our benchmark studies showed that HybridMine produced biologically accurate results for hybrid genomes compared to other well-established software. As proof of principle, we carried out a comprehensive structural and functional annotation of complex yeast hybrids to enable system biology prediction studies. HybridMine is developed in Python, Perl, and Bash programming languages and is available in GitHub.